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Abstract 

, We review the Weyl- Wigner formulation of quantum mechanics in 

VjQ ' phase space. We discuss the concept of Narcowich- Wigner spectrum and 

! use it to state necessary and sufficient conditions for a phase space func- 

tion to be a Wigner distribution. Based on this formalism we analize 
the modifications introduced by the presence of boundaries. Finally, we 
discuss the concept of environment-induced decoherence in the context of 
■ the Weyl- Wigner approach. 
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1 The Weyl- Wigner formulation of quantum me- 



^ chanics 

The Wiener function was introduced in the 1930's [T]. Wi gner was trying to 
derive quantum corrections to Boltzmann's kinetic equation. With this aim in 
■ mind he had to define some quantum counterpart to the Liouville density p c i. 

Let us recall some basic facts about this distribution. p c \ is assumed to be 
a real, normalized, non-negative, C°° function in phase-space. Here, we assume 
a flat 2n-dimensional phase-space T*M ~ R 2 ". The expectation value of some 
observable A £ C ,OD (T*M; R) is then given by: 

<A> N = J dx J dp A(x,p)p c i[x,p). (1) 

Moreover, the probability distribution that a measurement of A yield the value 
a £ R is evaluated according to: 

V (A = a) = { dx { dp 5 (A(x,p) - a) p cl (x,p). (2) 
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Finally, the time evolution of both (1) and (2) can be computed in the Heisen- 
berg or Schrodinger pictures: 

< A(t) >= J dx J dp A(x,p,t)p c i(x,p) — J dx J dp A(x,p)p e i(x,p,t) 

V (A(t) =a) = Jdxjdp S {A(x,p, t) - a) p el [x, p) = (3) 

= f dx f dp 6 (A(x, p) - a) p ct (x,p, t) 

where A(x,p, t) and p c i(x,p,t) are the solutions of Newton's and Liouville's 
equations, respectively: A — {A, H} p , p c i — {H,p c i} p . 

Wigner's purpose was to establish some correspondence rule such that the 
expectation value of some hermitian operator A be given by: 

< A >= J dx J dp A(x,p)F(x,p), (4) 

in analogy with classical statistical mechanics. The objects A and F would be 
the phase-space counterparts of A and of the density matrix p, respectively. It 
is important to stress that the choice of F and A is not unique due to different 
possible choices of operator orderings. In this paper we shall adopt Wigner's 
choice, namely Weyl's rule pP-Jlj. We define a C-linear map - the Weyl map W 
- which attributes to an element A of the quantum algebra of linear operators 
Aq an element A in the "classical" algebra Ac according to: 



W(A) = (^\ J dtj dr, Tr{A(q,p)e~ l ^-">p}e l t x +"iP. 



(5) 



This map appears naturally as a replacement (q,p) < — > (x,p), once the operator 
has been written as a linear combination of completely symmetrized polynomials 
of q and p. 

An important consequence of this map is that it introduces a noncommu- 
tative product (the starproduct 0]) and a bracket (the Moyal bracket 0) in 
phase space via the relations: 



W 



(A-fy=A*B, =^[A,B] M , (6) 



where A = W(A) and B = W(B). From (5,6) one can obtain the explicit 
expressions: 

A* B = Ae^^B = A-B + 0(h) 

t -\ (?) 
[A, B] M = lA sin B = {A, B} p + O (ti 2 ) 

where J is the " Poissori'' operator: J = [■§-■§ — 7r~ir ) > the derivatives d and 



dq dp dp dq J ' 

d acting on A and B, respectively. From (7) we realize that the starproduct 
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and the Moyal bracket are ^-deformations of the usual commutative product 
of functions and of the Poisson bracket, respectively. The noncommutativity 
of the starproduct, on the other hand, is manifest, if we consider the relation 
A* h B = B*- h A. 

It is also interesting to remark the following property: 

J dx J dp A(x,p) * B(x,p) = J dx J dp A(x,p) ■ B(x,p). (8) 

Thus far, we established the relation between the observable A and the phase 
space symbol A — W(A). Similarly, we have F(x,p) — (2^K)~ n W (p). This 
object is called the Wigner function and is the quantum mechanical counterpart 
of the classical Liouville measure p c ;. If the system is in a pure state with 
p = \tp >< ip\, then the Wigner function reads: 

F(x,p) = j^y; J dr, e-^* (x - r,h/2) i>{x + V h/2) . (9) 

It is worth summarizing some of the properties of Wigner functions. They are 
real, smooth and normalized as the classical Liouville measures p c i. Likewise, 
expectation values of operators are computed according to the formula (4), 
reminescent of the classical expression (1). However, and contrary to p c i, the 
Wigner function takes on negative values |2J as is illustrated by the following 
simple example. Let us consider, the simple harmonic oscillator, with Weyl 
symbol: 

H(x,p) = ^- + ^mLu 2 x 2 . (10) 
The Wigner function of the first excited state, then reads: 

1 f 4H(x,p) \ ( 2H(x,p) \ 

We conclude that F\{x,p) is negative inside the ellipse: H(x,p) < 

The fact that Wigner distributions are not, in general, positive defined is not 
surprising. One would otherwise be entitled to interpret them as true probability 
measures in phase-space, which in turn would violate Heisenberg's uncertainty 
principle. This is why quantum distributions in phase space are called " quasi- 
distributions" . In any case, the marginal distributions are bona fide probability 
distributions. From (9) one gets: 

dp F(x,p) = \^(x)\ 2 > 0, J dx F(x,p) = \4>(p)\ 2 > 0, (12) 

where ip(p) 1S the Fourier transform of ip(x). Moreover, if Fx, F2 are the Wigner 
functions associated with the states \ip >, \(f> >, respectively, then one gets: 



dx I dp Fx(x,p)F 2 {x,p) = < V# > r > 0. (13) 
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The previous inequality remains valid for any pair of Wigner functions F± , F2 
whether pure or mixed, and it is the basis of the probabilistic interpretation in 
quantum phase space. In particular, if F a (x,p) is the Wigner function associated 
with the eigenstate \ip a > of some hermitian operator A with eigenvalue a, then 
the probability distribution that a measurement of A on a system in a state 
F(x,p) yield the eigenvalue a is given by 1 : 

V (A = a) = (2itK) n J dx J dp F{x,p)F a {x,p). (14) 

It is interesting to remark that, since the projector \tp a >< ip a \ satisfies A\ijj a >< 
V*o| = \tpa >< ipa\A = o-\ip a >< ipa\, then upon application of the Weyl map, 
one obtains |10j: 

A(x,p) * F a (x,p) = F a (x,p) * A(x,p) = aF a (x,p). (15) 

One calls this a stargenvalue equation and F a is said to be a stargenfunction of 
A(x,p) = W(A) with eigenvalue a. This equation admits the formal solution 

M 

(2nh) n F a (x,p) = 5, (A(x, P )~a) = -L / dk t £W**)-) ) (16) 

where ef = 1 + B + B * B/2 + B * B * B / (3!) + • • • is the non-commutative 
exponential. The distribution <5* is called the non-commutative delta function 
because of its similarities with the ordinary Dirac delta function. From (14) and 
(16), one gets: 

V (A = a) = J dx J dp F(x,p)5* {A(x,p) - a) , (17) 

which is the quantum mechanical analog of (2). 

We now turn to the dynamics of the system. The density matrix is a solution 



of the von Neumann equation: ihdp/dt = 



H,p 



If we multiply this equation 



by (2tvK) " and apply the Weyl map, we obtain the Moyal equation: 

F(x,p,t) = [H(x,p),F(x,p,t)] M = {H(x,p),F(x,p 1 t)} p +0(h 2 ) (18) 
which is a deformation of the classical Liouville equation. 

In this work, we will apply the previous formalism to systems with boundaries 
(section 3) and to decohering systems interacting with a heat bath (section 4). 
In the latter situation we shall derive some known results using the concept of 
Narcowich- Wigner spectrum (which is discussed in section 2) and suggest some 
further developments. 

1 We assume for simplicity a nondegenerate, continuous spectrum. More general cases are 
addressed exhaustively in ref. 1141 



4 



2 Narcowich-Wigner Spectrum 



So far, starting from ordinary quantum mechanics, we derived the necessary set 
of rules to make physical predictions in quantum phase space. But if Wigner 
quantum mechanics is to be a self-consistent autonomous formulation of quan- 
tum mechanics in its own right, we should be able to reverse the arguments, 
i.e. starting from Wigner quantum mechanics we should be able to re-derive 
Schrodinger quantum mechanics. 

The obvious questions to pose are: (i) given some smooth, real and normal- 
ized phase space function F(x,p) how do we know whether it is a bona fide phase 
space representative of a quantum state? (ii) if F(x,p) represents a quantum 
state, how do we know whether it is a pure or a mixed state? 

To answer the first question, we start by stressing that not all functions 
F(x,p) satisfying the aforementionned set of properties are Wigner functions. 
Take for instance 

F(x,p) = ^e- aH ^\ (19) 

Z7T 

where H (x, p) is given by (10) and a is some constant with dimensions [energy) 1 . 
This function is C°°, normalized and real. However, the product of dispersions, 
A22; • A2P = violates Heisenberg's uncertainty principle whenever a > 
in which case it cannot represent a quantum state. 

To state the set of necessary and sufficient conditions > we first need to 
define the notion of pure state in quantum phase space. We say that a real, 
smooth, normalizable function F{x,p) is a pure state whenever it satisfies the 
differential equation: 

— lnZ(x,y) = (~) _lnZ(x, y ), (20) 

where Z(x,y) — J dp e lyp F(x,p). Alternatively, the previous equation can be 
replaced by the more algebraic condition: 

F(x,p)*F(x,p) = ^F(x,p). (21) 

From the previous equation we obtain, 

dx J dp F 2 (x,p) = JL (22) 

We thus construct the set T pure of all pure states in quantum phase space. The 
space of all states T is the set of real, normalizable, smooth functions F(x,p) 
such that: 

J dx J dp F(x,p)F pure (x,p) > 0, VF pure € T pur e- (23) 

States in phase space satisfy the inequality: 

J dx J dp F 2 (x,p) < JL, (24) 
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the equality holding only for pure states (cf.(22)). Mixed states also admit 
an algebraic representation, of which (21) is the particular expression for pure 
states. They are positive elements of the algebra ^3], i- e - F(x,p) belongs to T 
iff there exists a smooth, square integrable function g{x,p) such that 2 : 



F{x,p) = g(x,p) *g(x,p). (25) 

The set of necessary and sufficient conditions for F(x,p) to be a Wigner function 
stated above can be rephrased in other terms more suited for certain applica- 
tions. But first we introduce the notion of symplectic Fourier transform of some 
function f(x,p): 

/(a) = [dz f{z)e^ a \ (26) 



where we used the compact notation z T — (x,p), a T = (u, v) and cr(z, a) is the 
symplectic form: 

a(z, a) — z T J a = (x,p) ( ^ J ~\ f " ~\ — xv — up. (27) 
The inverse symplectic Fourier transform is: 

s{z) = W^S dz /» eMa ' z) - ( 28 ) 

A well known property of Fourier transforms is that they turn convolutions, 

(f*g)(z)= [ dz' f(z-z')g(z') 7 (29) 



into products: 

(Ffg)(a) = f(a).g(a). (30) 
We then state the following [23] 

Definition 1: Let /(a) be a continuous function on the dual of the phase space. 
Then / will be termed of 77-positive type if, for every m € N and any set of m 
points {ai, • • • , a m } in the dual of the phase space, the m x m matrix M with 
entries 

M jk = f( aj -a k )e ir >^ ak ' a ^ 2 , (31) 
is self-adjoint and non-negative. 

We are now in a position to state the set of necessary and sufficient contions for 
F to belong to T equivalent to the conditions (23) |22 | -|26 | . 

2 The function g(x, p) is defined up to ^multiplication by an "unitary" phase space function, 
i.e. g(x,p) and u(x,p) *g(x,p) yield the same state provided u(x,p) *u(x,p) = 1 
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Definition 2 (KLM Conditions): A function F(z) is a Wigner function iff its 

symplectic Fourier transform F(a) satisfies the KLM (Kastler, Loupias, Miracle- 
Sole) conditions: 

(i) F(0) = 1, 

(ii) F(a) is continuous and of ?i-positive type. 

Similar conditions also characterize classical states, i.e. the Liouville measures. 
We just need to replace ^-positive by 0-positive type, which is equivalent to 
F(z) being nonnegative everywhere in phase space (Bochner's theorem). 

More generally, we shall denote by Narcowich- Wigner (NW) spectrum W (F) 
of the function F(z) the set of real values n such that F is of ^-positive type. 
In particular, if F is a Wigner function, then Ti £ W(F). Another important 
point is that if a function F(z) is of 77(7^ 0)-positive type, F(0) = 1 and F{a) is 
continuous, then F{z) is continuous, square integrable and vanishes at infinity. 

For j] > we can define an 77-starproduct: A B = Ae lr,J / 2 B. It then 
follows that if F is of ^-positive type, then there exists some smooth square 
integrable function g such that: F{z) — g(z) * n g(z). Given the properties of 

the starproduct, we have: F(z) = g(z)*^ v g(z) — (^g(z)j *^ v g(z). We conclude 

that if 77 £ W (F), then also —77 £ W (F). Moreover, it can be shown that (for 
?7 ^0): 

dz F 2 (z) < -Lr. (32) 
27r|ry| 

In summary, the best way to think about rj(> 0)-positive functions is to regard 
them as Wigner functions where h is replaced by r\. 

Now, suppose that F(z) is a pure state Wigner function and that 77 £ W (F). 
From (22) and (32) we have: J dz F 2 (z) = ^7- < 2^f\- We have thus proved 
the following: 

Lemma: If F £ T pure then W (F) C [-h, h]. 

The previous lemma concerning the NW spectrum of a pure state can be further 
refined: 

Theorem 1: Let F £ T pure . If F is a Gaussian then W (F) ~ [— fr, h] , otherwise 

W(F) = {~h,h}. 

The proof of this theorem will be given elsewhere |16|. As an example, Narcowich 
computed the NW spectrum of the Gaussian, 



det A 



= \/ ex P [- ( z ~ z of A ( z ~ z o)\ > (33) 

where A is a real, symmetric, positive defined 2n x In matrix. It was shown 
that F is of Ty-positive type iff the matrix B = A^ 1 + in J is non- negative. 
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In particular, if A = cr x ld (in suitable units), where c is some real positive 
constant, then W (F) = [— c, c]. 

Another interesting property of NW spectra stems from the convolution of 
distributions. In fact, let the Schur product of two matrices A and B be the 
matrix M with entries Mjk = AjkBjk- Then it can be shown that if A and B 
are non-negative, then so is M. From the definition of NW spectrum, we can 
prove the following |25j : 

Theorem 2: Let W(F),W (G) and W (F * G) be the NW spectra of F, G and 
of the convolution F*G. Then we have: W (F) + W (G) C W (F-kG). 

In the previous relation, the set on the left-hand side is W (F) + W (G) = 
{»7 + xl»?e W(F), x e W(G)}. 

We now consider a simple application of the concept of the NW spectrum. 
Let Fq(z) be some real, continuous function with non-empty NW spectrum 
W(-Fb). What are the set of conditions that Fo should satisfy so that its con- 
volution with any Wigner distribution F G F yield another Wigner function? 
Using the properties of NW spectra it is easy to derive a set of sufficient condi- 
tions. Whether they are also necessary remains unproved. The conditions are 
as follows |25| : 

(i) Let F be of 0-positive type. If F € J 7 , then ±h € W (F). From theorem 
2, we conclude that {— h, ft,} C W (Fq * F). Consequently, Fq ★ F is a Wigner 
function for any F G T . 

(ii) Let ±2ft G W (F)) and let F G T . Similarly, from theorem 2, we obtain 
{±ft, ±3ft} C VV (F * F). Again, F * F is a Wigner function for any F G F. 

These results have two interesting applications. First, suppose that Fq is also 
a Wigner function (±h G W (F))). We conclude that Fq * F is of 0-positive 
type for any F G F. This means that under conditions (i) and (ii), Fq -k F is a 
non-negative Wigner function. This has an interesting physical interpretation. 
If Fo is taken to be a state associated with an apparatus, then Fo * F can be 
regarded as being the state that results when one uses that apparatus to make 
a measurement on a system in the state F [251 ■ The second application will be 
postponned to section 4. 

3 Wigner Functions with Boundaries 

Consider a particle of mass m in the infinite potential well |15j , confined to the 
interval a < x < b. For simplicity we shall assume Dirichlet boundary condi- 
tions. In ordinary quantum mechanics we solve the free Schrodinger equation 
and impose the boundary conditions thereafter. The normalized solutions are 



8 



(a = -b = -L/2): 



ip n (x) = y/2/Lsm[(mr/L)(x + L/2)], (34) 

inside the box and n — 1, 2, 3, • • • The energy eigenvalues are E n = h 2 ir 2 n 2 / (2mL 2 ) . 
If we compute the associated Wigner function using (9), we obtain for —L/2 < 
x < L/2: 

F n (x,p) = <^ cos (2^) • sin [f (| - 1*1)] + 

sin [wl (P L + (I ~ N)] + (35) 

sin [it (p £ ~ (t ~ N)] 

It is easy to verify that this expression satisfies Dirichlet boundary conditions 
F n {—L/2,p) = F n (L/2,p) — 0, Vp S fR. However, and contrary to the 
wavefunction, it also obeys Neumann boundary conditions dF n /dx (—L/2,p) — 
dF n /dx(L/2,p) = and also d 2 F n /dx 2 {-L/2,p) = d 2 F n /dx 2 (L/2,p) = 0, Vp G 
IR. But the most important aspect of this problem is the following. We expect 
the free particle in phase space to obey the stargenvalue equations: 

2 2 

— * F n (x,p) = F n (x,p) * — = E n F n (x,p), (36) 

with the same eigenvalue E n . However, a straightforward calculation shows that 
F n (35) is not a solution of (36). This may be linked to the non-local nature 
of the ^-product, which entails that the presence of the boundaries is felt well 
inside the bulk. Also notice that since e l ^^ h implements spatial translations and 
translation invariance is broken when there is a boundary, we conclude that the 
Heisenberg algebra is no longer valid. Remember that the Heisenberg algebra 
is at the heart of the Weyl quantization. 

So the problem we want to solve is the following 3 : 

(i) How do we modify eq.(36) to accommodate the boundary effects? 

(ii) How do we modify the Moyal equation for the dynamics, when a boundary 
is present? 

(iii) How do we choose the boundary conditions satisfied by the Wigner function? 

First of all let us go back to the definition of the Wigner function (9). In the 
confined case, ij)*{x — y)ip{x + y) vanishes outside of the parallelogram a < 
x — y < b, a < x + y < b. If Xq = (a + b)/2, then F(x,p) can be defined 
sectionwise, for a < x < xq: 

F(x,p) = F 1 (x,p) = l~r a dy e- 2l Py' h r{x - y)1>(x + y), (37) 
nri J a - X 

3 For an alternative recent approach we refer to 1201 . 
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and for xq < X < b: 

r-b 



F(x,p) = F 2 (x,p) = -1 / dy e- 2i ™/ h r(x ~ yW* + v). (38) 

Since the wavefunction or its derivatives are possibly discontinuous at the bound- 
aries and to avoid possible misinterpretations, the previous integrals are de- 
fined as improper: ff_ ° stands for lim c _ >a + f x _ ° and, likewise, J f stands for 

lim c ^ b - f°_ ■ From definitions (37,38) it is easy to check that F is continuous 
at x — xq and satisfies Dirichlet boundary conditions irrespective of the bound- 
ary conditions satisfied by the wavefunction. The boundary conditions on -0 
do play a role in the derivatives of F, namely: dF/dx(a + ,p) = 2|-0(a + )| 2 /(7r7i) 
and dF/dx(b~ ,p) = — 2\ip(b~)\ 2 / (irh). And so, if tp obeys Dirichlet boundary 
conditions, then: 

dF, , , dF n OF _ dF , , . 

to ( ° ' P) = to r ' P) = °» to {X ° ' P) = to (X ° (39) 



but also: 



d 2 F, . . d 2 F „ . d 2 ^, _ , <9 2 F 



fe2 (« -P) = fe^'P) = °> to (Xo "' P) = to (X ° (40) 

However, 0(a+,p) = %W {a+)\ 2 / {nh) £ 0. 

We are now in a position to compute the boundary corrections to the star- 
genvalue equations f36l|T5): 

Theorem 3: If the wavefunction ip is a solution of the Shrodinger equation 

- h 2 /(2m)ij/'(x) + V(x)ip(x) = E^{x), (41) 

with Dirichlet boundary conditions, then the associated Wigner function F(x,p) 
solves 4 : 



(42) 



+ L S '( X ~ b ) * F ( x " CP) = EF{x,p). 
Moreover F{x,p) obeys the boundary conditions: 



j dp F(a + e,p) = J dp F(b- e,p) = 0. (43) 



Proof: A straightforward calculation shows that (for a < x < xq): 

^ + V(x)\*F 1 (x,p)=EF 1 (x,p)+B(x,p) 1 (44) 



4 The infinitesimal parameter e is to be taken to + once the starproducts and the integra- 
tions have been evaluated. 
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where the boundary contribution is: 

B ( X >P) = -2L e ' 2w(a ' X)/h {nf^( 2x - a~ma+) + 
+tp'*(2x - a-)ip(a+) + tP*{2x - a -)tp'{a+)} . 
For Dirichlet boundary conditions, the previous equation reads: 



(45) 



B D (x,p) = --^ e - 2tp{a - x)/h ip*{2x - a-)V'(a+). (46) 
z7rm 

A simple calculation shows that this can be expressed in terms of the Wigner 
function as: 



B D (x,p) = lim — ( dk e tk{x - a) kF 1 {x + e,p~Tik/2). 
e^o+ Airm J 



(47) 



The parameter e appears in the previous expression because we want the deriva- 
tive ip'(a + ) and not Tp'{a~) (remember that ij)' is discontinuous at a). If we 
expand B D in powers of h, we get: 

h 2 

B D (x,p) = lim 5'(x -a)* F(x + e,p). (48) 

2m e—o+ 

Performing a similar calculation for Fi in xq < x < b and assembling all the 
results we obtain (42). That F obeys the set of boundary conditions (43) is a 
trivial consequence of the relation lim e _ + J dp F (x ± e,p) = (x ± ) | 2 . 

Some remarks are now in order. First of all, it can be checked that eq.(35) 
is a solution of the boundary stargenvalue equation (42). Secondly, we can 
perform Baker's converse construction. This means the following. Suppose that 
some real, normalized function F is solution a of cq.(42) and that it satisfies 
the boundary conditions (43). Then, we can show that there exists a square 
integrable function ip (unique up to a phase) which (i) obeys Dirichlet boundary 
conditions at a and b, (ii) is a solution of the Schrodinger equation (41) and (iii) 
is related to F(x,p) via the Weyl transform (9). Thirdly, a similar calculation 
for the dynamics yields the following boundary Moyal equation: 

%(x,p,t)=\£; + V(x),F(x,p,t)\ + 

(49) 

+£ i S '( x - b )> F ( x - t'P^U - & t - o),F(x + e,p,t)] M , 
as expected. 

4 Decoherence 

The paradigm of environment induced decoherence consists of a Brownian par- 
ticle interacting with a heat bath. The Brownian particle may be regarded as 
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an open system, the behaviour of which will usually be described (under general 
conditions) by master equations with a local time. In a typical situation |18j 
of a dust particle interacting with air molecules or radiation, the dust particle 
can then be described by a Markowian master equation for the reduced density 
matrix 5 : 

S-5SS^]-f' , ^-«l- (50) 

Such an equation appears frequently from an interaction with the environment 
in cases where friction is negligible. The strength of the coupling is given by the 
parameter D. The first term in (50) would lead to unitary spreading, whereas 
the second term would lead to nonunitary localisation. The time scale for de- 
coherence is: 

rah 

to = \J^- (51) 

If we change to the Weyl-Wigner representation of quantum mechanics, we 
conclude that the corresponding Wigner function satisfies the Fokker-Planck 
equation: 

dF _ p dF D d 2 F 

~dt ~~m~fa +1 2~dp Y ' ( ' 

If decoherence occurs, the general expectation is that the negative parts of the 
Wigner function will gradually be smoothed out in the course of time. Many 
examples support this expectation |18j . 

In ref. p2| Diosi and Kiefer proved a much stronger statement, namely that 
the Wigner function becomes strictly positive after a certain finite decoherence 
time to- This is very different from the behaviour of the density matrix whose 
non-diagonal elements become zero only asymptotically. This is important be- 
cause the positivity of the Wigner function is usually regarded as a necessary 
condition for classicality. 

Another relevant aspect of this mechanism is the fact that the time tjj , after 
which the Wigner function becomes strictly positive, is independent of the initial 
Wigner distribution. 

The results of that paper were later generalized in ref. |19j , where a Markovian 
master equation, 



dp if i \ 



dt h 



H,p 



2h ^ 



fapLt-L^Ljp-pLlLj), (53) 



was considered with a set of Lindblad operators {A?} linear in the Brownian 
particle's position and momentum and representing the interaction with the 
environment. The conclusions were identical: the Wigner function becomes 
strictly positive after a finite time irrespective of the initial condition. 

Let us rederive the results of [TB] using the formalism of section 2. The 
Green's function for the system is the smooth function (for any t > 0) K t (z\z r ), 

5 This is the operator which results when the environment's degrees of freedom have been 
traced out from the closed system's density matrix. 
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solution of the Fokker-Planck equation (52) with initial condition Kq (z\z r ) = 
8(z — z'). We can then write: 

F(z,t) = J dz' K t (z\z')F (z'), (54) 

where F n (z) — F(z,0) is the initial Wigner distribution. A straightforward 
calculation leads to the following solution. The propagator K t {z\z') is of the 

form (33), with n = 1, z$ = (x' + and: 

6m 2 _ 3m \ 

(55) 

_ 3m __2_ ] 
IF Dt I 

Since the propagator is of the form K t (z\z') — Q t (z — Zo), we may, alternatively, 
rewrite (54) in the form: 

F(z,t) = (Q t *F )(z_ t ), (56) 

where Z- t = [x + j^,p)- From (33), we know that 77 beleongs to the Wigner 
spectrum of fi t (z) iff the matrix, 

/ _2Dt? 

B t =A; 1 +i V J=l , (57) 

V ^ + iv -2Dt J 

is non-negative. We conclude that the Wigner spectrum of Slt(z) is: 

Dt 2 Dt 2 ' 



W(fi t ) = 



m\/3' m\/3. 

It follows that n t (z) becomes a Wigner function once t>to, where: 



(58) 



4 3m 2 h 2 

The Wigner function ^(2, t) (for t > t^j) is then the convolution of a 0-positive 
Wigner function fit(z) and another Wigner function F (z) (cf(56)). From the 
analysis in section 2, this means that F(z,t) is non-negative from tu onwards. 
Eq.(59) is the result obtained by Diosi and Kicfer. Notice that to is independent 
of the initial Wigner distribution Fq(z). 



How general is this mechanism? The reason why tu is independent of F (z) 
stems from the fact that it is the propagator that enters the space T of Wigner 
distributions at t = tn and this only depends on the dynamics of the system, 
i.e. on the master equation. Thus, suppose that the Wigner function at time 
t > is given by eq.(56), where zt is some trajectory such that zq = z and fit (2) 
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is a propagator. Assuming that the set of conditions derived by Narcowich 
and analyzed in section 2 are both necessary and sufficient, we may define 
the set of propagators K. as follows. We say that fit(z) belongs to IC if (i) 
fit (2) is smooth for t > 0, (ii) fio(z) = 5(z), (iii) fi t satisfies the semi-group 
rule (fi t *fi t ,)0) = fi t+t '0), for t,t' > and, finally, (iv) S W(fi t ), or 
±2?i € W(fit). The latter conditions ensure that the convolution of fit with 
an initial Wigner distribution Fq(z) yield another Wigner distribution. The 
propagator is the solution of some master equation with initial condition (ii). 
Suppose that at to the propagator becomes a Wigner function, i.e. h g W (fit) 
for t > to. Then, concomitantly the Wigner function F(z,t) becomes strictly 
positive. From this construction it is clear that tjj cannot depend upon the 
initial distribution Fq(z). 

What we propose is the following. Instead of trying to derive kinetic equations 
for open Markowian systems, which is notoriously difficult for non-quadratic 
interactions (the only ones considered in refs.JS|> we ma y instead look 

for propagators in the set K, which, at a given time, enter the set T . We may 
even consider situations where the propagator enters T at time tc and leaves 
it at a later time tQ. This is presumably what happens in ref.|17|. This is the 
programme of a future work [lfi| . 
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